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Abstract 

In this paper, we study decoherence in continuous-time quantum walks 
(CTQWs) on one-dimension regular networks. For this purpose, we assume 
that every node is represented by a quantum dot continuously monitored by 
an individual point contact(Gurvitz's model). This measuring process in- 
duces decoherence. We focus on small rates of decoherence and then obtain 
the mixing time bound of the CTQWs on one-dimension regular network 
which its distance parameter is Z > 2. Our results show that the mixing 
time is inversely proportional to rate of decoherence which is in agreement 
with the mentioned results for cycles in pU [37]. Also, the same result is 
provided in [38] for long-range interacting cycles. Moreover, we find that 
this quantity is independent of distance parameter 1(1 > 2) and that the 
small values of decoherence make short the mixing time on these networks. 
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1 Introduction 



Quantum walk(QW) as a generalization of random walk(RW) is attracting 
great attention in many research areas, ranging from solid-state physics [lj 
to quantum computing [2]. Experimental implementations for the quan- 
tum walks have been presented in [31 01 [5]. In recent years, two types 
of the quantum walks exist in the literature: the continuous-time quantum 
walks(CTQWs) [0 El El El [TO] and the discrete-time quantum walks(DTQWs) 
[IH CE3 El Q31 ESI US]. The relationship between the CTQWs and the 
DTQWs has been considered in H3 HB1 US] . The CTQWs have been stud- 
ied on star graph [20], [21] , on direct product of cayley graphs [22] , on quotient 
graphs [23] , on odd graphs [23] , on trees [25] and on ultrametric spaces [26] . 
All of these articles have focused on the coherent CTQWs. The effect of 
decoherence in the CTQWs has been studied on hypercube [27j [28], on 
cycle [29j . on line [301 131] , on iV-cycle [32] and on long-range interaction 
cycles |38| . Here, we study the CTQWs on one-dimension(lD) networks 
with diatance parameter I > 2 which can be constructed as follows [33]: we 
construct an one dimensional ring lattice of N nodes, each node of which 
is connected to its 21 nearest neighbors(7 on either side). The structure 
of one-dimension regular network with N = 8 and I = 3 is illustrated in 
Fig. 1. One-dimension regular networks have broad applications in various 
coupled systems, for example, Josephson junction arrays [53] . small-world 
networks [35j and synchronization [36J. In our paper, the network nodes are 
represented by identical tunnel-coupled quantum dots(QDs). The walks are 
performed by an electron initially placed in one of the dots. An individual 
ballistic one-dimension point-contact is placed near each dot as "detector" 
which its resistance is very sensitive to the electrostatic field generated by 
electron occupying the measured quantum dot. Decoherence is induced by 
continuous monitoring of each network node with nearby point contact. We 
focus on small rates of decoherence, then calculate the probability distribu- 
tion and the mixing time bound of the CTQWs on one-dimension regular 
network with distance parameter I > 2. Our analytical results show that 
small decoherence can make short the mixing time in the CTQWs. The 
same result was produced for cycles in [29] and for long-range interacting 
cycles in |38| . Moreover, we show that for small rates of decoherence, the 
mixing time is independent of distance parameter 1(1 > 2). 
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Figure 1: One-dimension regular network with N = 8 and I = 3 



This paper is organized as follows: 
In Sec. 2, we briefly review the properties of CTQWs on one-dimension 
regular networks. In Sec. 3, we study the decoherent CTQWs on the under- 
lying network. We assume that the rate of decoherence is small and obtain 
the probability distribution, analytically in Sec. 4. The bound of the mixing 
time and its physical interpretation are provided in Sec. 5. Conclusions and 
discussions are given in the last part, Sec. 6. 

2 CTQWs on ID regular network 

The properties of network is well characterized by the spectrum of adjacency 
matrix of associated graph. The network adjacency matrix(^4) is defined in 
the following way: A{j = 1 if nodes i and j are connected and otherwise 
Aij = 0. The Laplacian is defined as L = A — D, where D is a diagonal 
matrix and Dj j is the degree of vertex j. Classically, the continuous-time 
random walks(CTRWs) are described by the master equation [TJ [39] 



where p^j (t) is the conditional probability to find the walker at time t and 
node k when starting at node j. The transfer matrix of the walk, T, is 
related to the adjacency matrix by T = —jL. For the sake of simplicity, we 
assume that the transmission rate 7 of all bonds to be equal. The formal 
solution of Eq. (1) is 



The quantum-mechanical extension of the CTRW is called the continuous- 
time quantum walk(CTQW). The CTQW is obtained by replacing the Hamil- 
tonian of system with the classical transfer operator, H = —T [1, 39, 40J. 
The Hamiltonian matrix H for one-dimension regular network is written as 
the following form [33] 




(1) 



p k>j (t) = (k\e It \j). 



(2) 
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—2m, if i = j; 

Hij = (i\H\j) = { +1, if i = j±z,z e [1,1]; (3) 
0, Otherwise, 

that the basis vectors \ j) associated with the nodes j span the whole acces- 
sible Hilbert space. In these basis, the Schrodinger equation(SE) is 



ij t \j)=H\j), (4) 

where we set m = 1 and H = 1. The Hamiltonian acting on the state \j) can 
be written as 



H\j) = -(21 + l)\j) + \j + z ), z e Offers (5) 

which is the discrete version of the Hamiltonian for a free particle moving on 
a lattice. It is well known in solid state physics that the solutions of the SE 
for a particle moving freely in a regular potential are Bloch functions |41|, 142] . 
Thus, the time independent SE is given by 

H\^ n )=E n \^ n ), (6) 

where eigenstates \<& n ) are Bloch states. The periodic boundary conditions 
require that $ n (iV) = <J? n (0), where <& n {x) = {x\^ n ). This restricts the 9- 
values to 6 = where n = 0, 1, . . . , N — 1. The Bloch state |$ n ) can be 
expressed as a linear combination of the states \j) localized at nodes j, 



N-l 

I 



1 N ~ X 



n v^ 

Substituting Eqs. (5) and (7) into Eq. (6) , we obtain the eigenvalues of the 
system as 

E n = -2l + 2Y,cos{je n ). (8) 
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The time evolution of state \j) starting at time to is given by U(t,to)\j), 
where U(t,to) = exp(—iH(t — to)) is the quantum mechanical time evolution 
operator. Hence, the transition amplitude akj(t) from state |j) at time 
to state \k) at time t is 

a k j(t) = (k\e' im \j). (9) 

Applying E n and |3> n ) to represent the nth eigenvalue and eigenvector of H, 
the classical and quantum transition probabilities between two nodes can be 
written as 

PuAt) = Y. e ~ tEn ^)^\i)' ( 10 ) 

n 

n ktj (t) = \a k>j (t)\ 2 = \Y,e- UE "(k\$ n )(<S> n \j)\ 2 . (11) 

n 

3 The Decoherent CTQWs on ID regular network 

In this section, we investigate decoherence induced by the point contact(PC) 
detector measuring the occupation of one of the quantum dots(QDs) in a 
double-dot system. The measurement process is shown schematically in 
Fig. 2. We assume all electrons to be spin- less fermions and the tunneling 
between the PCs and the QDs to be negligible, but we take into account 
Coulomb interaction between electrons in the QD and the PC. We start with 
writing the Hamiltonian for the entire system. The total Hamiltonian is 

N-l 

H = H s + ^2 (Hpcj + Hi n tj), (12) 

3=0 

where H s , H pCt j and Hintj would be identified next. 

Note that in this paper, quantum walk is defined over an undirected graph 
with iV nodes that each node is labeled by an integer n £ [0, iV— 1]. Also, we 
assume that the quantum walker has no internal state (i.e. simple quantum 
walker), so that we can describe its dynamics by a Hamiltonian of form |43j : 

h s = Yl A a + - Yl ^(*)4<v> 

(13) 

ij 3 
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Figure 2: Fig. 2(a) shows point contact detector j monitoring the electron 
in dot j and Fig. 2(b) shows point contact detector j when electron is placed 
in dot j + 1. [ii j and chemical potentials of left and right reservoirs 

of jth point contact, respectively. Ej is the on-site node energy. 

where cj(cj) is creation (annihilation) operator and the walker at node 

j corresponds to the quantum state \j) = ct 1 0} . The first term in Eq. (13) 
is a 'hopping' term with amplitude Ajj(i) along the edge ij between nodes 
i and j; the second term describes 'on-site' node energies Ej(t). We assume 
the hopping amplitude between connected sites to be constant and drop on- 
site energy terms (ie., Ej = 0, Vj). Also, for convenience, we renormalize 
the time so that it becomes dimensionless [H]. Thus, the simple quantum 
walker Hamiltonian has the form: 

N-l I 

^ = iEE(4+^'+4 c ^)- ( 14 ) 

j=0 z=l 

The tunneling Hamiltonian H pc describing this system can be written as 
h pcJ = E + E ^i a r,A,j + E n lr,j(4,j & r,j + al, 3 ai,j), (15) 

l r l y r 

where a\ j(aij) and dj j{a r j) are creation(annihilation) operators in the left 
and right reservoirs of point contact j, respectively. Also, Eij and E r j are 
the energy levels in the left and right reservoirs of detector, and £li r j is the 
hopping amplitude between the states E[j and E r j. We assume that the 
hopping amplitude of jth point contact is S7/ r ,j when an electron occupies 
the left dot, and it is Qi r j when an electron occupies the right dot. Hence, 
we can represent the interaction term as 

Hint = ^2 ^irjc'jCjial jarj + o^ajj). (16) 

l,r 
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where (5£li r j = 0/ r j — fijrj). For simplicity, we assume that the hoping 
amplitudes are weakly dependent on states Eij and E r j, so that £li r ,j = ^, 
SCli r j = 5Q and Wj(jU r ,i) = Pi(Pr)- Gurvitz in [55] applied the Bloch- 
type equations for a description of the entire system with the large bias 
voltage(/Uj — /i r ). He showed that the appearance of decoherence leads to 
the collapse of the density matrix into the statistical mixture in the course 
of the measurement processes. Using Eq. (5), this analysis for our model 
results in the following equation for the reduced density matrix 

I 

TkPj,k{t) = |[ (^'> fc + 2 ~~ Pj+z,k)] ~ r (! ~ S j,k)Pj,k 
z=-l 

(17) 

i 

= j[^2(Pj,k+z - Pj+z,k + Pj,k-z - Pj-z,k)\ - r(l - Sj : k)Pj,k- 

2 = 1 

4 Small Decoherence 

In this section, we assume that the decoherence rate T is small as (TiV <1) 
and consider its effect in the CTQWs on one-dimension regular network. 
For this end, we make use of the perturbation theory of linear operators ,as 
mentioned in [29], and rewrite Eq. (17) as the perturbed linear operator 
equation 

N-l 

Jt P«At) = E (* x (5) + tf£S)W*)> ( 18 ) 

fi,u=0 

where a,(3,[A,v run from to N — 1. Also, L^ V X and U^ v 2 which are 
the row (/i, u) and column (a, (3) elements of N 2 x iV 2 matrices of L and U 
respectively, are defined as 

1 1 

L (a!f3) = j[ E i 8 oL,^p,v-z ~ 8*,n-zfy,v)], (19) 



= -r<W<W(i - <W). (20) 

Also, for our case, the initial conditions are 
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(21) 



Now, we want to obtain the eigenvalues and the eigenvectors of L + U. 
For this aim, we study the perturbed eigenvalue equation [29] 

(L + U)(V + V) = (X + X)(V + V), (22) 

that V is the corresponding eigenvector with eigenvalue A i.e. LV = XV. 
Applying first-order perturbation theory of quantum mechanics, one can get 

A = V ] UV. (23) 

We assume that e\ is the eigenspace with eigenvalue A and some of eigen- 
vectors of L (i.e. {Vk '■ k E /}) span it. For the uniform linear combination, 
we can obtain A as following 

X = Y J V]uV k . (24) 
fee/ 

The solution of Eq. (18) is obtained by dropping the terms Vj [29] 

P(t) = £ et(a+X) £c^i- (25) 
For unperturbed linear operator L, we have 

N-l 
fi,u=0 

that A (m>n ) is 

v— v . .TTz(n + m), .irz(m — n). , 
\m,n) = E Sm (^x7 ) Sm ( N ) ( 27 ) 

2 = 1 

and V M 18 



S 



Using Eq. (24), one can obtain 

U im ,n),(m>,n>) = {V^^UV^' ^ 



W X^ 1 ~ 6a ^ ex P(-^-[( m ' - m ) a + ( n ' - n ) 6 D (29) 

(a,b) 

Y5 m i m 5 n i n J^^[(m'— m)+(n' ~~n)](modN),0 

In what follows, we first find the degenerate eigenvalues of Eq. (27) and 
then calculate the eigenvalue perturbation terms: 

(a) Diagonal element (m = n): 

Since U is diagonal over the corresponding eigenvectors, there is not such 
degeneracy in our case |29j . For these eigenvalues, the correction terms are 



given by Eq. (29): \ m . m) = -T^fi. 

(b) Zero(m + n = 0(modN)): 

Since the corresponding eigenvectors can not display in the linear combina- 
tion of the initial state p(0), this degeneracy is irrelevant to our problem |29j . 

(c) Off-diagonal elements: 

By Eq. (29), the off-diagonal elements are non-zero if 
m + n = m! + n'(modN). 

To find degenerate eigenvalues satisfying the relation(c), we make divide the 
problem into two separate states as follows: 

The state 1 = 1: This state is equal to a cycle network, for which Eq. 
(27) reduces to \ m ^ n ) = sin( n ( n + m > ) s i n ( m±p) ) _ \( m ^ = \ m ,^ results 

■ /Tr(m+n)\ , ■ /TT<m'+n')\ i ■ /Trim— n) \ , ■ /7r(m'— n') \ mi 

m sm(— — '-) = ± sm( v N — '-) and sm(— y — — '-) = ± sm( v N — '-). inus, 
we have 

' m = n > + N/2, n = m' + N/2, for N/2 < m < N - 1, N/2 < n < N — 1; 
m = n > - N/2, n = m! - N/2, for < m < N/2, < n< N/2; 
m = n' + N/2, n = ml - N/2, for N/2 <m<N - 1,0 < n< N/2; 
m = ri - N/2, n = m' + N/2, for < m< N/2, N/2 < n < N - 1. 
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The correction terms to these eigenvalues are A( m n ) = — r ^ N N 2 ^ ■ 

The state I > 2: In this case, the structure of cycle network can be de- 
stroyed by 21 additional bonds in the network. 

\ra,n) = \m',n') implies to 

E 7rz(n + m) irz(m — n) ^-^ /Kz(n' + m') .irzim' — n 1 ), 
sm( )sm( ) = ^sin( )sm( . 

2=1 2=1 



One can see that there is not any degeneracy for this state. Thus, cor- 

•y-r 

N 



rection terms to these eigenvalues are \ m n = —T ^ N ^ 



The mixing time bound for the state I = 1 was provided in [29} 137]. 
In the following, we focus on the state I > 2 (one-dimension regular network 
under condition I > 2) and in the end, compare our results with [29\ 137] 's 
results for cycle. 

Based on the above analysis, Eq. (25) can be written as 

p(t) = ^(<W,o + S m+n>N )V^ (30) 

(m,n) 

and from Eq. (21), we have 

= E ^ (m ' n) - (31) 



N 

(m,n) 



Thus, the full solution is 



(m,n) 



(32) 



x exp[^(ma + n(5)]. 

The probability distribution Pj(t) of the quantum walk is specified by 
the diagonal elements of the reduced density matrix, i.e. 

Pjit) = 7f + W Yl 0- ~ $[m+n]{modN),o) X [e~ F ^*] 
(m,n) 

I (33) 

r ,TTz(m + n) .TTz(m — n).. . . . 

x exp [it 2^ sm( — )sm( — )] x exp[^(m + n)j\. 

2 = 1 
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5 Mixing time 



There are two distinct notions of mixing time for quantum walks in the 
literature: 

Instantaneous mixing time: Instantaneous mixing time is defined as the first 
time instant at which the probability distribution of the walker's position 
is e-close to the uniform distribution [46J. Thus, the instantaneous mixing 
time is 

M inst)t = min{t : \\Pj(t) - —\\ tv < e}, (34) 

where here we use the total variation distance to measure the distance be- 
tween two distributions P, Q: \\P — Q\\t v = Yli 1-^ (*) ~~ Q(i)\- 
Now by Eq. (34), we calculate an upper bound on the Mi nst:t 

, „ . . -, . pS-i,. i v-^ / r. v— ^ . ,-Kzim + n) . ,irz(m — n) 

\PjV)-7t\ =e~ r — t \w ^(exp^^sm( 1 - , )sm( V - ; )] 

(m,n) z=l 

xexp[^(m + n)j]) - 



By noting that 

^-s irz(m + n) irz(m — n) .2m '■ . 
exp^2^ sm ( m ) sm ( )J x exp[— (m + n)jj < 1, 

2 = 1 



one can obtain 



Thus, we have 



m)-^\<e-^\ (36) 



N ~ 1 1 

^|^-(t)-^|<iVe- r V*. (37) 

3=0 

Based on the above definition, the upper bound of instantaneous mixing 
time can be obtained in the following way: 



Ne~ r — * < e , (38) 
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and therefore, 



IN 1 

M inst , e <-ln(-)[l + ¥ - I ]. (39) 

The Mi ns t te for cycles was provided in [29], that is 



IN 2 
M inst , e <- ln(-)[l + ¥ — (40) 

These relations show that the instantaneous mixing time for one-dimension 
regular network with distance parameter 1(1 > 2) is shorter than the one for 
cycle network. Also, this quantity is independent of distance parameter I . 

Average mixing time: To define the notion of average mixing time of CTQWs, 
we use the time-averaged probability distribution, i.e. Pj(t) = - Jq Pj(t)dt. 
The average mixing time measures the number of time steps required for 
the time-averaged probability distribution to be e-close to the limiting dis- 
tribution 1471. i.e. 



M ave , t = min{t\\/T > T : \\Pj(t) - —\\ tv < e}. ( 1.1 ) 



In the following, we want to obtain the lower bound of average mixing time 
for one-dimension regular network(Z > 2). 
Applying Eq. (36) for large N ^> 1, we have 

IIP (t) _ III < 1 1 f T ( e ~rt + ]_ )dt _ 1 1 = J_r x _ -rri (42) 
ii n > N w _ i T v t N ) N \ rr i j \ ) 

Summing over j to calculate the total variation distance, we have 

^(l-e- rT )<6. (43) 

Then we assume that TT > 1 (since N » 1 and TN < 1, thus T > N) 
and achieve 

N 

Mave,e > ~ ; (44) 

re 

which is similar to the average mixing time bound produced for cycle in [37] . 
Physical interpretation: 

As mentioned in Sec. 3, we assumed the hopping amplitude between all 
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of connected sites to be equal. Also, we supposed that every point-contact 
detector which is coupled to states localized on the corresponding node, 
measures the position of the particle in space of graph. Since these detectors 
identify the path of the walker took, quantum interference which is the result 
of an uncertainty in the path, is then lost and therefore the mixing time is 
independent of distance parameter I [48]. However, one notes that this 
explanation might not be broadly accepted. 

6 Conclusion 

We studied the effect of small decoherence on one-dimension ring lattice of 
N nodes in which every node is linked to its 21 nearest neighbors (I > 2 
on either side). In our investigation, this network was represented by the 
system of identical tunnel-coupled quantum dots. As the detector, we used 
the point contact in close proximity to one of the dots. For description of 
the entire system, we applied the Bloch-type equations. Then, we calculated 
the probability distribution and the mixing time bound. We showed that 
the mixing time bound is independent of parameter 1(1 > 2). Also, we 
observed that this quantity is inversely proportional to rate of decoherence, 
as mentioned in [29} I37j . Hence, decoherence can make short the mixing 
time on these networks. Moreover, we found that the instantaneous mixing 
time for one-dimension regular network is smaller than the one for cycle 
network. 
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